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Abstract 

Mutations in the human DNA methyl transferase 3A {DNMT3A) gene are recurrently identified in several hematologic 
malignancies such as Philadelphia chromosome-negative myeloproliferative neoplasms (MPN), myelodysplastic syndromes 
(MDS), IVIPN/IVIDS overlap syndromes and acute myeloid leukemia (AML). They have been shown to confer worse prognosis 
in some of these entities. Notably, about 2/3 of these mutations are missense mutations in codon R882 of the gene. We 
aimed at the development and validation of a novel easily applicable in routine practice method for quantitative detection 
of the DNMT3A p.R882C/H/R/S mutations bead-based suspension assay. Initial testing on plasmid constructs showed 
excellent performance of BNA(NC)-modified probes with an optimal hybridization temperature of 66°C. The method 
appeared to be quantitative and showed sensitivity of 2.5% for different mutant alleles, making it significantly superior to 
direct sequencing. The assay was further validated on plasmid standards at different ratios between wild type and mutant 
alleles and on clinical samples from 120 patients with known or suspected myeloid malignancies. This is the first report on 
the quantitative detection of DNMT3A R882 mutations using bead-based suspension assay with BNA(NC)-modified probes. 
Our data showed that it could be successfully implemented in the diagnostic work-up for patients with myeloid 
malignancies, as it is rapid, easy and reliable in terms of specificity and sensitivity. 
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Introduction 

Somatic mutations in the DNA metliyltransferase 3A 
(DNMT3A) gene have been reported as recurrently associated 
with myeloid malignancies. The frequency of mutations varies 
between different entities. In adult acute myeloid leukemia (AML) 
DNMT3A mutations are found in 14-34% of cases from different 
series [1], 5-15% of MDS cases [2], 10% of chronic myelomono- 
cytic leukemia (CMML) patients [3], 5.7% of primary myelofi- 
brosis (PMF) patients [4], 12% of cases with systemic mastocytosis 
[5] and about 18% of T cell acute lymphoblastic leukemia (T- 
ALL) cases [6]. Approximately 2/3 of the DMMT3A are missense 
mutations at codon R882 [7]. Most of the current literature 
supports the notion that the presence oi DMMT3A mutations is an 
adverse prognosis biomarker at least in adult AML and may be an 
important parameter of the integral molecular genetics profiling in 
these patients [8,9,10]. These lines of evidence necessitate the 
development of methods for reliable detection of DJVMT3A 
mutations that are easily applicable and affordable for routine 
use together with the testing for other mutations. The most 
frequently reported techniques used for DJVMT3A mutations 
include direct sequencing, high-resolution melting analysis and 



next generation sequencing. Sanger sequencing is particularly 
well-established technique [7,11] for the identification of previ- 
ously unreported mutations but its relatively low sensitivity (~10- 
20% mutant allele burden) may be problematic for detection of 
low frequency somatic mutations. High resolution melting (HRM) 
technique has also been adapted for detection of DNMT3A 
mutations and is also good for screening for unknown mutations in 
a single tube format with a sensitivity of about 4% [12,13]. 
However, it requires weU-estabhshed standards and eventually 
verification of the result through direct sequencing. Targeted 
amplicon resequencing on next generation sequencing (NGS) 
platforms has also been used for DMMT3A mutations [14,15]. The 
obvious advantages of NGS as the digital allele burden output, 
theoretically very high sensitivity, the possibility for identification 
of novel mutations are currently limited by the costly equipment 
and the necessity for a strong bioinformatic support. Simple 
approaches such as restriction fragment length polymorphism 
(RFLP) based assay have also been developed [16], whose 
relatively low sensitivity could be overcome by couphiig with 
quantitative PGR (qPCR)[17]. Another less frequently used 
technique for verification of DMMT3A mutations confirmation 
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Table 1. Primary sequence of the tested oligonucleotides and the predicted melting temperatures (Tm). 



Modification 


Probe designation/Codon cKiange 


Sequence 


Tm 


LNA 


m (CGC) 


5 '-CGCCAAGCGGCTCAT-3 ' 


68°C 




R882C (CGC->TGC) 


5 '-CGCCAAGCAGCTCAT-S ' 


66 C 




R882H (CGC->CAC) 


5 '-CGCCAAGTGGCTCAT-S ' 


66°C 




R882P (CGC->CCC) 


5'-CGCCAAGGGGCTCAT-3' 


67"C 




R882S (CGC->AGC) 


5 ' -CGCC A AGCTGCTC AT-3 ' 


66°C 


BNA{NC) 


WT (CGC) 


5 '-CGCCA+AGCG+GCTCATGrT-B ' 


70 C 




R882C (CGC->TGC) 


5 '-CGCC+A+AGCAG+CTC+ATGIT-B ' 


70°C 




R882H (CGC->CAC) 


5 '-CGCC+A+AGrG+GCTC+ATGIT-S ' 


70 C 




R882P (CGC->CCC) 


5 ■ -CGCC A+AGGGG+CTCATGTT-3 ' 


70°C 




R882S (CGC->AGC) 


5 ■ -CGCC-l-A-l-AGCTG-l-CTC+ATGTT-3 ' 


70 C 



The character "+N" denotes the position of any BNA(NC) modified nucleotide. The positions of LNA modification are proprietary property of Exiqon, Denmark and are, 

therefore, not disclosed. 

doi:10.1371/journa!.pone.0099769.t001 



was denaturing high-performance hquid chromatography 
(DHPLC) [18]. 

In our view a particularly suitable for routine diagnostics could 
be the detection of somatic mutations in hematologic malignancies 
using a microsphere-based suspension array. It allows single tube 
quantitative detection of mutations in one or several amplicons in 
a mid-throughput format, making it a rapid and affordable 
method. Here we report the development and validation of a 
sensitive multiplexed bead-based suspension assay for quantitative 
detection of DMMT3A R882 mutations. The method appeared to 
be applicable to clinical samples from patients with myeloid 
malignancies. 

Materials and Methods 

Ethics statement 

The study was conducted in accordance with the principles of 
the Declaration of Helsinki. Written informed consent was 
obtained from all patients. Blood and bone marrow sampling as 
well as molecular testing were part of the routine diagnostic 
procedures approved by the Institutional Review Board (IRB) at 
Alexandrovska University Hospital, Sofia, Bulgaria. 

Patients 

A total of 120 peripheral blood or bone marrow samples of 
consecutive patients with known or suspected myeloid malignan- 
cies were collected between January 2010 and June 2013. Patients 
were classified according to the WHO criteria as follows: Acute 
myeloid leukemia (AML) (n = 21). Chronic myeloid leukemia 
(CML) (n= 1), Myelodysplastic syndrome (MDS) (n = 6), Myelo- 
proliferative neoplasm (MPN) (n = 81). Overlap MPN/MDS 
(n — 5), and others with suspected but unproven MPN (n = 6). 

DNA extraction 

AH samples were collected using sodium citrate-containing 
blood sampling tubes (BD Biosciences) and stored at room 
temperature for no more than 4 hours before processing. Genomic 
DNA was extracted from whole blood using a iPrep (Invitrogen) 
automated system. Genomic DNA samples were stored at — 20°C 
before further analyses. 



DNMT3A exon 23 sequencing 

Exon 23 of the human DMMT3A gene was amplified from 
genomic DNA as reported before [5] using the following primers: 
DNMT3A (forward): 5'-TCCTGCTGTGTGGTTAGACG and 
DNMT3A (reverse): 5'-ACAGAAAACCCCTCTGAAAAG. Am- 
plification in 25 |J,1 reactions included the following: 100 ng 
genomic DNA; 1.5 U Taq polymerase (Invitrogen); 3 mM MgCl2; 
0.2 mM dNTP mixture; and 10 pmol of each primer. The 
amplification conditions were as follows: 95°C for 5 min; 35 cycles 
of 95°C for 30 sec, 58°C for 40 sec, 72°C for 2 min and final 
extension at 72°C for 10 min. PGR products were purified using 
Exo-SAP (Applied Biosystems, USA) and sequenced bidirection- 
ally with the same set of primers. The sequencing reaction was 
performed in a final volume of 1 0 |J,L using 2 |J,L of the purified 
PGR product, 3.2 pmol of one the primers and 2 |J,L of Big Dye 
terminator cycle-sequencing kit v3.1 (Apphed Biosystems, USA), 
The sequencing program was 25 cycles of denaturation at 96°G 
for 10 s, annealing at 50°G for 5 s and extension at 60°C for 
4 min. The sequence detection was conducted using the ABI 
Prism 3100 Genetic Analyzer (Applied Biosystems, USA). 
Sequences alignment and analysis were performed using Se- 
quencher v. 5.0. (GeneGodes, USA). 

Bead-based assay 

The bead-based assay was performed as described before 
[19,20,21]. Briefly, the exon 23 DMMT3A fragment was amplified 
from either genomic or plasmid DNA samples using a 5'- 
botinylated forward primer. The same primers and PGR 
conditions as described above for the sequencing analysis were 
apphed. Genotyping was performed by direct hybridization with 5 
LNA- or BNA(NC)-modified oligonucleotide probes, specific for 
the wild type or the mutant alleles. LNA-modified oligonucleotides 
were designed and synthesized by Exiqon, Denmark, and the 
BNA(NG)-modified ohgonucleotides were designed by us and 
synthesized by BioSynthesis, USA. The sequences and the melting 
temperatures (Tm) of all oligonucleotides are shown in Table 1 . AU 
probes were synthesized with 5 '-amino group and 20 bases as a 
spacer sequence for the purpose of covalent binding to the 
carboxylated microspheres (Luminex, USA). The coupling of the 
amine modified probe to the carboxylated surface of the beads was 
performed using a standard carbodiimine-coupling procedure 
[22]. Mixtures of 5 sets of coupled microspheres were prepared by 
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Figure 1. Fitted linear standard curves from tKie test assays for determination of the best performing probe set and the optimal 
hybridization temperature. 

doi:1 0.1 371 /journal.pone.0099769.g001 



combining equal volumes of each set of the microspheres. 
Approximately, 160 beads of each type/|J,l were used for the 
analysis of one sample. Five (0.1 of PGR product; 20 [il 
Hybridization buffer (Wakunaga Pharmaceuticals, Japan); 3 [ll 
Bead mixture and 2 |J,1 of streptavidin-phycoerythrine (SAPE) 
(W akunaga Pharmaceuticals, Japan) were combined per well of a 
96-well Thermowell plate and hybridization was performed at 
62°C, 64°C, 66°C or 68°C for 30 min in a thermocycler (Applied 
Biosystems, USA). After incubation, 75 |J,1 Washing buffer 
(Wakunaga Pharmaceuticals, Japan) were added to each well. 
The supernatants were removed after plate centrifugation at 
4000 rpm for 2 min. Finally, the microspheres were resuspended 
in 75 |J,1 Washing buffer (Wakunaga Pharmaceuticals, Japan) and 
acquired on a LabScan^"" flow platform (Luminex, USA). A 
minimum of 100 events per bead region of interest was collected. 
For each set of reactions a background control of a sample 
containing only the respective sets of microspheres was included. 
The background median fluorescence intesity (MFI) values were 
subtracted from the MFI values for each sample. The resulting 
values were used for the calculation of the relative fluorescence 
indices for each mutant allele according to the following equation: 
Index (mutant allele;) = MFI (mutant allele,)/ [MFI (mutant 
allele/) + ... +MFI (mutant allele ^) + MFI (wild type allele)]. 

Determination of the cut-off index values of the assay 

Sensitivity and specificity of the assay for each mutant as well as 
the optimal threshold values for mutant allele indices were 
determined using Receiver Operating Characteristic (ROC) 
analysis [23]. All indices for positive and negative samples for 
each mutant from five independent runs of the assay were used for 
the analysis. The optimal cut-off point was determined by the 



Fisher's test implemented through the R language-based web tool 
Cutoff Finder (http://molpath.charite.de/cutofr/) [24]. Analysis 
of the sensitivity in terms of percent mutant allele burden was 
performed by plotting the range of indices for each mutant at 
concentration 1%, 2.5% and 5% and 12.5% versus the already 
determined optimal cut-off value. 

Results 

Determination of the optimal probes set and 
hybridization conditions 

The development of any hybridization-based assay for detection 
of single nucleotide changes (e.g. SNPs, point mutations) is 
particularly challenging because of the need to design highly 
specific probes that are able to discriminate between the difiFerent 
alleles with high accuracy. The task becomes even more 
complicated when more than two alleles are to be discriminated 
in a multiplex fashion. The design of the probes is also influenced 
by the overall nucleotide content of the fragment of interest and in 
some instances even single nucleotide changes can lead to a high 
variability in the melting temperatures. So, the primary goal of the 
optimal probe set design is to select probes with high and uniform 
Tm. Traditionally, this could be achieved through variation of 
probe length, which is a time and labor-consuming process. To 
overcome this limitation a number of chemically modified bases 
have been used to design oligonucleotide probes. For the purpose 
of developing a bead-based assay for detection oi DMMT3A R882 
mutations we tested the performance of two sets of oligonucleotide 
probes with synthetically modified nucleotides. The first set 
comprised of f generation bridged nucleic acids (BNA), known 
also as locked nucleic acids (LJSrA)-oligonucleotides, i.e. they 
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contain one or more 2'-0,4'-C-methylene-P-D-ribofuranosyl 
monomer(s) [25]. These probes were designed and synthesized 
by the LNA patent liolder Exiqon, Denmark. The second set of in- 
house designed probes included third generation BNA ohgonu- 
cleotides, known also as generation bridged nucleic acids 
(BNA(NC)), i.e. they contain more than one 2'-0,4'-aminoethy- 
lene-P-D-ribofuranosyl monomers [26]. The design of the 
BNA(NC) ohgonucleotides was performed following the published 
recommendations (http:/ / www.biosyn.com/ faq/BNA-general- 
design-guidelines.aspx) and they were synthesized by BioSynthesis 
Inc., USA. 

As shown in Table 1 the BNA(NC) probes had uniform Tm of 
70°C, whereas the LNA probes had slighdy lower Tm values 
varying between 66°C and 68°C. Based on our previous 
experience with Luminex-based hybridization assays[21] we 
expected the optimal performance of each set of probes to be 
around 2°C-4°C below the expected Tm. Indeed, to determine 
the optimal hybridization conditions for each set of probes we used 
dilutions of 100%, 50%, 20%, and 0% of each mutant plasmid in 
the wild-type plasmid. For the LNA probes set we tested 62°C and 
64°C hybridization temperatures and for the BNA(NC) probes we 
tested 66°C and 68°C hybridization temperatures. Mutant allele 
indices for each sample were calculated and linear regression 
model was fitted for each mutant at every hybridization 
temperature (Fig. 1). The parameters for every linear model fit 
are presented at Table 2. The BNA(NG) probes set performed 
better than the LNA probes and the former was selected for 
further testing. The 66°C and 68°C hybridization temperature 
levels tested for the BNA(NC) probes were actually identical in 
terms of the parameters of the fitted linear models (Table 2). We 
selected the 66°C for all subsequent runs of the assay as at this 
level the raw MFI values were higher than those at the 68°C level. 

Determination of the index threshold value and assay 
sensitivity 

To determine the mutant allele index threshold values we used 
Receiver Operating Characteristic (ROC) analysis. ROC curves 
were based on the analysis of the index values from plasmid 
samples with wild type only (true negative samples) and wild type 
and mutant alleles samples (true positive samples) as described 
above. ROC curves for all mutants are shown at Fig. 2. The 
determined areas under the curve (AUC) were indicative for an 
excellent performance of the assay with all values being over 0.9 as 
follows: P.R882C (0.94), p.R882H (0.98), p. R882P (0.94), and p. 
R882S (0.97). The corresponding specificity and sensitivity at the 
selected cut-off index values are shown in Fig. 2. In order to 
determine the corresponding sensitivity of the assay for each 
mutant in terms of mutant allele burden we assayed plasmid 
standards with 1%, 2.5%, 5%, and 12.5% of each mutant plasmid 
on wild type background. The sensitivity was defined as the lowest 
mutant allele burden whose range of index values was above the 
determined respective threshold. As shown in Fig. 3 the 
determined sensitivity was 2.5% for all DMMT3A mutants. 

Concordance of the bead-based assay with direct 
sequencing using genomic DNA 

We further tested the performance of our bead-based assay on a 
set of genomic DNA samples from 120 patients with known or 
suspected myeloid malignancies. We identified a total of two R882 
missense mutations (one p.R882C and one p.R882Q) (Fig. 4). 
Both cases were AML. The results from the bead-based assay for 
all cases were compared with direct sequencing of exon 23 of 
DMMT3A. The sequencing data confirmed the presence of the two 
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Figure 2. ROC curves for the determination of tKie overall performance of the assay and the optimal cut-off value for each mutant's 
index. 

doi:1 0.1 371/journal.pone.0099769.g002 



R882 mutations but also identified another unreported heterozy- 
gous synonymous mutation p.Qf)86Q^ (c.26.'i8G>A). This patient 
was also diagnosed with AML. All other cases were negative for 
R882 mutations on both the microsphere-based assay and direct 
sequencing. This straightforward comparison validated the appli- 
cability of our assay for the detection of R882 mutations in clinical 
settings. 

Discussion 

Recent advances in genomic profiling technologies showed that 
genes encoding enzymes involved in epigenetic regulation are 



frequentiy affected by the mutational process in myeloid malig- 
nancies [27]. DMMT3A is a key gene for the epigenetic regulation 
of human cells as the encoded enzyme is known to be a major 
player in the de novo DNA methylation at CpG sites [28]. 
DNMT3A mutational status has recentiy been included in several 
studies on the prognostic power of integrated molecular profiling 
in AML [9,10,29]. Most of the studies used targeted sequencing of 
exon 23 of the gene spanning the mutational hotspot, codon R882. 
From a clinical standpoint the restriction of molecular testing for 
DNMT3A mutations only to the mutational hotspot, codon R882, 
is justified by three lines of evidence: i) approximately 2/3 of all 
reported DNMT3A mutations are affecting codon R882 (with 
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Figure 3. Determination of the sensitivity of tKie assay as expressed in percent mutant allele. The horizontal lines denote the threshold 

level for each mutant's index. 

doi:1 0.1 371 /journal.pone.0099769.g003 



P.R882H and p.R882C being the most frequent ones); ii) R882 
missense mutations perturb the normal molecular frinction of the 
enzyme; and iii) R882 mutations may differ in their prognostic 
effect from non-R882 mutations [18]. 

It was shown recently that the most frequent DJVMT3A 
mutation, R882H, disrupted the dimerization interface of the 
enzyme. This leads to a several fold drop down in its activity as a 
processive CpG methylase [30]. However, DM'IT3A R882 
mutations are not just hypomorphic and they seem to exert also 
dominant negative effect as the expression of the analogous mouse 
DNMT3A R878H mutations in embryonic stem (ES) cells 
containing endogenous DNMT3A or DNMT3B caused hypo- 
methylation [31]. Xu et al. [32] just showed that the transplan- 
tation of hematopoietic stem cells (HSCs) transduced retroviraUy 
with R882H mutant and transplanted to immunodeficient mice 
caused CMML disease with characteristic hypomethylation 
patterns. A number of retrospective clinical studies assessed the 
prognostic power of DMMT3A in AML. Only a few of them, 
however, compared R882 mutations with non-R882 mutations, 
but the conclusions regarding their prognostic power are not 
consistent. For instance, Gaidzik et al. [7] showed that R882- 
mutated cases may be prognostic for unfavorable OS vs. non- 
R882-mutated cases. Other studies, however, did not find such a 



difference between different types of DNMT3A mutations [33]. 
Taken together, the clinical and experimental evidence suggest the 
relevance of DMMT3A mutations detection in the work-up for 
AML (and eventually other related myeloid malignancies). 

To address the rapidly growing need for testing the DJVMT3A 
R882 in clinical settings we adapted the bead-based flow platform 
developed by the Luminex Corporation, USA. This system utilizes 
barcoding of various analytes recognizing molecules (antibodies, 
ligands, oligonucleotides) through coupling with fluorescently 
labeled microbeads [34]. In the cases of direct hybridization 
based method the amplified DNA fragment is either directly 
labeled with Cy3/PE or with biotin with subsequent streptavidin- 
phycoerythrine (SAFE) labeling. A number of recent reports 
demonstrated the applicability of this technique for detection and/ 
or quantitation of somatic mutations associated with malignant 
diseases such as KRAS and BRAF [35,36], MPMl, JAK2, MPL, 
IDHl/2 [19,20,21,37,38,39] and for detection of chromosomal 
translocations [40,41] and even for CpG methylation profiling 
[42] . The greatest advantage of this method is the opportunity for 
multiplexing in a single tube the detection of up to 100 analytes. A 
nice demonstration of that advantage for detection of cancer- 
associated mutations has recently been provided by Bando et al. 
[43] who performed multiplex single tube testing for 36 mutations 
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Figure 4. Comparison of the bead-based assay to direct sequencing in selected genomic DNA samples. The vertical arrows show the 
positions of the identified mutations. 
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in KRAS, BRAF and PIK3CA. Unfortunately, they did not 

provide any data on the analytical sensitivity of the assay. The 
greatest challenge in the development of such assays is the design 
of probes that could specifically discriminate between sequences 
differing in just one nucleotide base. In the last decade a working 
strategy to overcome this limitation appeared to be the use of 
synthetic nucleotides that can increase the specificity and the Tm 
of the oligonucleotide probes they are incorporated in. The LNA- 
containing oligonucleotides have been reported to be successfully 
coupled with Luminex beads for SNP or somatic mutations 
detection. Here, to the best of our knowledge we reported for the 
first time the use of BNA(NC)-containing oligonucleotides coupled 
with micro-beads and resolved on the Luminex platform. Our in- 
house designed BNA(NC) probe set appeared to be superior to an 
LNA probe set intended for the same purpose. The BNA(NC) 
incorporation allowed for more uniform Tm of the probes in the 
set resulting in better sensitivity and specificity as demonstrated by 
the very high linearity between the mutant allele; standards and the 
obtained mutant allele indices (Fig. 1 and Table 2). This 
BNA(NC)-modified probes superiority in comparison to the 
LNA-modified ones might not be considered surprising as this 
has already been shown for RNA silencing application [44] . 

We further used a state-of-the art technique to determine the 
index cut-ofiF value for presence of any of the specific mutant 
alleles. Notably, all the cut-off values translated in 2.5% mutant 
allele burden sensitivity, which is of clinical relevance, as this value 
would mean 5% heterozygous mutant carrier cells in the sample. 
In principle enhancement of the sensitivit)^ of the assay is possible 
through use of mutant-enrichment PGR protocols (COLD-PCR) 
or signal amplification techniques (e.g. use of PE-labeled anti- 
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